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ABSTRACT 

It has been shown recently that the dynamical I/-band mass-to- light ratios of compact 
stellar systems with masses from 10 6 M© to 10 s M Q are not consistent with the pre- 
dictions from simple stellar population (SSP) models. Top-heavy stellar initial mass 
functions (IMFs) in these so-called ultra compact dwarf galaxies (UCDs) offer an 
attractive explanation for this finding, the stellar remnants and retained stellar en- 
velopes providing the unseen mass. We therefore construct a model which quantifies 
by how much the IMFs of UCDs would have to deviate in the intermediate-mass and 
high-mass range from the canonical IMF in order to account for the enhanced M/Ly 
ratio of the UCDs. The deduced high-mass IMF in the UCDs depends on the age of 
the UCDs and the number of faint products of stellar evolution retained by them. 
Assuming that the IMF in the UCDs is a three-part power-law equal to the canonical 
IMF in the low-mass range and taking 20% as a plausible choice for the fraction of 
the remnants of high-mass stars retained by UCDs, the model suggests the exponent 
of the high-mass IMF to be ~ 1.6 if the UCDs are 13Gyr old (i.e. almost as old as 
the Universe) or rj 1.0 if the UCDs are 7Gyr old, in contrast to 2.3 for the Salpeter- 
Massey IMF. If the IMF was as top-heavy as suggested here, the stability of the UCDs 
might have been threatened by heavy mass loss induced by the radiation and evolu- 
tion of massive stars. The central densities of UCDs must have been in the range 
1O 6 -1O 7 M pc~ 3 when they formed with star formation rates of 10-100 M Q yr _1 . 
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1 INTRODUCTION 

Ultra compact dwarf galaxies (UCDs) are stellar systems in 
which 10 6 M Q to 10 s M Q of gas were converted into stars 
withi n a volume of some ten pc in diameter (IHilker et all 
19991: iDrinkwater et all l200d , 120031 : iPhillipps et all l200ll : 
Hasegan et al.ll2005| y~ If UCDs are ess entially the massive 



these properties indicate that UCDs once were among the 
most extreme star-forming regions in the universe. 

A fundamental function underlying star formation is the 
stellar initial mass function (IMF), £(m), 



dN oc £(m) dm, 



(1) 



end o f the globular clus ter sequence (|Mieske et all |2002| . 
120041 : iForbes et all 120081 ) , then this must have happened 
within a few Myr, so that the star formation rate would 
have been 10-1 00 M^yr -1 . Indee d , the enhancement in a- 
elements, that lEvstigneeva et all (|2007l ) found in most of 
the UCDs they examined, suggests a short time scale for 
the formation of their stellar populations. Taken together, 
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where dN is the number of stars with initial masses between 
m and m + dm. The IMF is the parent distribution for the 
mass functions of stars in star clusters |Kroupa fc Weidnerj 
2003) . These mas s functions are subject to statistical scatter 
jElmegreenl Il997l : iKroupal 120011 ) and have an upper mass 
limit determined by the mass of the gas cloud ou t of which 
the star cluster formed (|Weidner fc Kroupall2006l ). 

One of the most debated questions concerning the IMF 
is whether it is universal, i.e. independent on the conditions 
under which star formation takes place. This is not expected 
from a theo r etical point of view. lAdams fc Fatuzzol (|l996l ) 
and iLarsonl l|l998l ) suggest an increase of the characteris- 
tic masses of pre-stellar cloud cores with increasing ambient 
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temperature. iMurrav fc Lfnl (|l996h discuss interactions of 
pre-stellar clumps leading to mergers as a process in star 
formation. Their model predicts an increase of the mean 
stellar mass with the density of the star-forming region. At 
the transition from massive globular clusters (GCs) to UCDs 
(i.e. in the mass range between 10 6 and 10 7 M o ), encoun- 
ters between pre-stellar clumps must have been particularly 
important. Only about 100 times the diameter of the or- 
bit of Neptune is available for the mean distance between 
stars in the central parts of so me of these high-mass GCs 
or low-mass UCDs (see fig. 4 in iDabringhausen et al.ll2008l . 
hereafter DHK). If expansion due to mass- loss through gas 
expulsion and stellar evolution played a role during their 
youth, then the densities of UCDs would have been even 
higher at their birth. 

On the other hand, all observed resolved stellar popula- 
tions are consistent with having formed with the same IMF. 
This canonical IMF can be formulated as a two-part power 
law, 



with 

Qi = 1.3, 
a 2 = 2.3, 



f c (m) = kirn 



m 

0-1 < — < 0.5, 

~ M Q 

m 

0.5 ^ -— - ^ m mi 

Mo 



(2) 



where m max is a function of the natal stellar mass of an 
embedded star cluster at the ti me when star f ormat ion is 
over and £ c = for m > m max ljKroupall200l] . l200St ). The 
factors ki ensure that the IMF is continuous where the power 
changes. 

During the past years suggestions for the IMF not be- 
ing universal, but over-abundant in high-mass stars (top- 
heavy) under extreme conditions, have accumulated for dif- 
ferent types of stellar systems. These incl ude galaxies (e.g. 
Baugh et al"1l2005l . iNagashima et al.ll2005l and Ivan Dokkum 



2008]) . the Galactic bulg e and centre (e.g. [Ballero et a~ 2007 
amHManess et alll2007r) and Galactic globular clusters (e.g . 
ID'Antona fc Caloill2004l and IPrantzos fc Charbonne]|l200ah 
Especially Milky Way globular clusters (MWGCs) have 
been examined closely. Their stellar mass functions might 
have been altered strongly by early residual gas expulsion 
( Mar ks et al.l 120011) and ste l lar and dynamica l evolution 
Baumgardt fc Makinol 120031; iBorch et all 120071 ; iKruiissenl 
20081 ; Kruiissen fc LamersT 20081 ). but the observation of in- 



dividual stars in the MWGCs can still give clues on their 
IMFs; namely by interpreting the complex patterns of the 
element abund ances in MWGC st ars (e.g. the Na-O anti- 
correlation, see lGratton et aL 2004 for a review on the com- 
position of MWGC stars). These peculiarities are usually 
taken as evidence for self-enrichment, meaning that the last 
stars that formed in a particular MWGC contain material 
that has been processed by stars that formed earlier in the 
same cluster. 

Different theories on how exactly the process of self- 
enrichment took place have been brought forward: the 
metal-enrichment in subsequent stellar generations could be 
caused by the ejecta of m assive asymptotic-g i ant-b ranch 
stars, as sugges t ed e. g. by iDAntona fc CaloT l|2004h and 
ID'Antona et al.l l|2007h , or by the winds from very massive 
stars, as suggested e.g. by IPrantzos fc Charbonnell (|2006l ) 



and iDecressin et al.l l|2007h . Yet both approaches require a 
top-heavy IMF, although residual gas expul sion from mass- 
segre gated clusters alleviates this need (Dccrcssin et alj 
|200S| ). 

It was shown e.g. in DHK and lForbes et all (120081 ) that 
GCs and UCDs do not constitute two clearly distinguishable 
populations, if a sample that covers the whole mass interval 
from GCs to massive UCDs is considered. This suggests a 
close relation between GCs and UCDs. It therefore seems 
well possible that the peculiarities in the element abun- 
dances that are found for stars in massive MWGCs could 
as well be present in the even more massive UCDs. But the 
only nearby objects that may be considered as UCDs and 
can (like the MWGCs) be resolved into individual stars are 
u) Cen and (at least to some extent) Gl in M31. Such obser- 
vations indeed show the stellar content of these most massive 
star clusters (or low- mass UCDs) to have a spread of metal- 
licit ies and ages fe.g iMeylan et al.ll200ll . iKavser et aill2006l 
and lVillanova et al.ll2007l ). 

However, there is an alternative way to set constrains on 
the IMFs of the UCDs, namely by the comparison with sim- 
ple stellar population (SSP) models. Various authors thereby 
found that the UCDs tend to have higher dynamical U-band 
mass-to- light (M/Lv) ratios than expected for any possi- 
ble stellar population that formed with the canonical IMF 
jHasegan et al.ll2005l;lHilker et al.ll2007l ; iReikuba et"aT]|2007l ; 
DHK; lMieske et al.ll2008l ). 

This resul t could indicate the p r esence of non-baryonic 
dark matter (lHasegan et al. |_ | 2005l ; iBaumgardt fc Mieskd 
120081 ). However, Murray ( 20081 ) argues that both numeri- 
cal simulations and observations of dwarf spheriodal galaxies 
hint to dark matter densities that are far too low to influence 
the dynamics of UCDs. This strengthens the notion that 
the high M/Lv ratios of the UCDs are t he consequence of 
an IM F different from the canonical one. iMieske fc Kroupal 
(2008) discuss an over-abundance of low-mass stars (i.e. 
stars with high M/Lv ratios) as a possible cause for the high 
M/Lv ratios of the UC Ds. They make testable p redictions 



based on the CO-index (|Kroupa fc Gilmordfl994D . Comple- 
mentary to their approach, this contribution is dedicated to 
top-heavy IMFs as an explanation for the high M/Lv ratio 
of the UCDs, which is in this scenario the consequence of 
a large number of remnants from burnt-out stars in them. 
The possible need for a top-heavy IMF also in the context 
of the element anti-correlations in massive GCs, as outlined 
above, makes this approach particularly attractive. 

This paper is organised as follows. In Section[2j the data 
sample used in this work is introduced. Section describes 
the model that is constructed for the stellar populations in 
UCDs. The results suggested by this model for the IMF 
of intermediate-mass and high-mass stars are presented in 
Section [4] Some implications of these results are discussed 
in Section [5] We summarise and conclude in Section [6] 



2 THE DATA SAMPLE 

The present paper is based on the data of GCs and UCDs 
compiled in IMieske et al.l l|2008l ). their table 5, because the 
chosen sample fulfils two requirements necessary for what is 
done in the present paper: 
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(i) Estimates for the dynamical mass have to be available 
for the objects. 

(ii) Estimates of the global metallicity of the objects have 
to be possible. 

This sample is currently the largest and most updated sam- 
ple of its kind. We note however that the results in DHK 
are qualitatively unchanged, although the present sample 
has been revised and enhanced compared to the sample they 
use. 

The term 'dynamical mass' refers to a mass estimate 
that is based on the velocity dispersion of the stars in the 
stellar system (derived from spectral lin e widths) and the 
spatial structure of the stellar system (see lHilker et al 1 120071 
for details). The mass estimates are therefore independent 
from the observed total luminosities of the stellar systems. 

The metallicities of the stellar systems are of impor- 
tance for the present paper because of their influence on the 
luminosity of stellar populations. Knowing them is therefore 
essential for creating models of stellar populations with a 
certain M/Lv ratio, which is the focus of the present paper. 

Besides newly estimated quantities, table 5 in 
iMieske etT al. (2008) also comprises numbers that are taken 
from the previous literature, as documented in their paper 
for the masses but not for the metallicities. Details on the 
origin of the metallicity estimates for objects with masses 
^ 2 x 10 6 Mr are given in Tab. [2] of our paper. When 
IMieske et al. I ((20081 ) make their own metallicity estimate 
from the (V — I) colours of the stellar systems they use 
the relation 

[Fe/H] = 3.27(V - I) - 4.50 (3) 

(eq. 4 in iKissler-Patig et al.| [l998). This has been done for 
all objects in their sample with masses < 2 x 10 6 Mq, unless 
the stellar syst ems are MWG Cs for which the metallicities 
are taken from lHarrisI l| 19961 ) (private communication with 
S. Mieske). 

Following Mie ske et al.l l|2008l ). we take an estimated 
mass of ^ 2 x 10 6 Mq as an easy-to-handle criterion to 
categorise a compact stellar system as a UCD instead of a 
GC. This mass marks quite well the transition from objects 
with GC-like proper t ies to objects with UCD-like proper- 
ties ijMarastonl [20051 : IMieske et ail l200Sl : DHK), including 
the on average distinctively higher M/Lv ratios of the more 
massive objects. Note that the two-body relaxation time ex- 
ceeds a Hubble time fo r systems larger than 2 x 10 6 Mq 
jMieske fc Krou pa 2008), which has been proposed as the 
defining prope rty to distinguish galaxies from star clusters 
l|Kroupalll998l : DHK). 



3 A MODEL FOR THE STELLAR 
POPULATIONS OF THE UCDS 

We now construct a model for the stellar populations of 
the UCDs under the assumption that the deviations of their 
M/Lv ratios from the theoretical expectation for the M/Lv 
ratio of a stellar population with the canonical IMF are 
caused by an IMF that varies for intermediate-mass and 
high-mass stars. The actual shape of the IMF in the UCDs 
cannot be specified from resolved stellar populations so far. 
The purpose of the following can therefore only be to give an 
idea by how much the IMF must deviate from the canonical 



IMF in order to account for the mismatch between obser- 
vations and theoretical expectations for the M/Lv ratio of 
the UCDs. 



3.1 The model ingredients 

The problem of modelling a stellar population with a M/Lv 
ratio equal to an observed value can be formulated as 

Ml 

Lr 



T v = 0, 



(4) 



where M m is the total mass of the model population, L m is 
its luminosity in the V-band and Ty is the observed M/Lv 
ratio of a stellar system. M m and L m depend on various 
parameters, such as the assumed age of the population, the 
shape of its IMF and the chosen model for stellar and cluster 
evolution. L m additionally depends on the metallicity. These 
dependencies will be formulated below, along with the as- 
sumptions that are made for the model presented here. 



3.1.1 The IMF 

The IMFs of the UCDs are connected to their present-day 
mass functions in the simplest way possible, because of their 
median two-body relaxation times, t r h, which are of the or - 
der of a Hubble time or larger (DHK; IMieske et aT1l2008l ). 
The timescale on which a stellar system dissolves depends 
on the tidal field strength, but can be expected to be many 
t r h, so that the stellar populations of UCDs are practically 
unaltered by dynamical evolution. This stands in contrast 
to GCs, whose are much shorter and therefore can have 
experienced significant dynamical evolution since their for- 
mation (also see Section |4. 2. 2[) . 

We introduce a family of IMFs for the model stellar 
populations of the UCDs: 



with 

Qi = 1.3, 
q 2 = 2.3, 

a :i £ R, 



£ p i(m) = kiva 



0.1 < — <0.5, 

°- 5 ^l£ <1 ' 
m 

1 ^ ^ ^max , 

Mq 



(5) 



where m max is the upper mass limit for stars. These IMFs 
will be referred to as the 'three-part power-law IMFs'. They 
are equal to the canonical IMF except for their slope above 
1 Mq. We assume that the UCDs have formed with a three- 
part power-law IMF. 

Upper mass limits of 100 Mq and 150 Mq are consid- 
ered. The upper mass limit of 100 Mq equals the upper mass 
limit assumed in the simple stellar population (SSP) models 
which are used in this paper (see Section [3. 1.2[) . These are 
the same stellar population models DHK took as a refer- 
ence when they found that the M/Lv ratios of a significant 
majority of the UCDs tends to be higher than model pre- 
dictions for the canonical IMF. This mass limit is however 
not in agreement with the uppe r mass limit for stars in ver y 
massive star clusters given by IWeidner fc Kroupal l|2004h . 
lOev fc Clarke] l|2005h and iFigerl l|2005h . which is close to 
150 Mq . Therefore this more realistic upper mass limit is 
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considered as well. It turns out that the results are affected 
surprisingly little by the upper mass limit of the IMF (Figs. [2] 
and [4] below) . 

Note that the lower mass limit of the IMF neglects the 
exist ence of brown dwarfs. T his is probably unproblematic, 
since iThies fc Kroupal (|2007l) showed that a combined mass 
function of brown dwarfs and stars shows a discontinuity. In 
the case that the low-mass IMFs of the UCDs are compa- 
rable to the ones in Galactic open star clusters (as assumed 
to be the case here), their results suggest that brown dwarfs 
contribute only a few percent to the total mass of the UCDs. 

The formulation of the IMF in UCDs given in eq. [5] at- 
tributes a possibly enhanced M/Lv ratio of a UCD solely 
to a top-heavy IMF, i.e. to a large population of stellar rem- 
nants that would have to be expected in such a case. Note 
however that the assumption of a bottom-heavy, Salpeter- 
Massey-like IMF in UCDs is currently an equally valid ap- 
proach to explain their M/Lv ratios (cf. figure 10 in DHK). 
Observations to test the hypothesis of a bottom-heavy IMF 
in UC Ds using a method proposed in iMieske fc Kroupal 
|2008f ) are underway. 



3.1.2 Simple stellar population models 

A simple stellar population (SSP) is defined as a popula- 
tion of stars of the same age and metallicity. Various au- 
thors have set up grids o f mo dels of such popu lations, e.g. 
iBruzual fc Chariot! <|2003h and IMarastonl (|2005l ). The alter- 
ation of the stellar mass function due to dynamical evolu- 
tion is not considered in these grids; only stellar evolution 
changes the mass spectrum of the stars in the model popu- 
lations. 

The most closely examined object in the sample of 
UCDs used here, uj Cen, is known to have several stel- 
lar sub-populations of differe nt ages and metallicities, i.e. 



uj Cen is not a SSP (e g.lmiker fe Richtlerll2000l ; iHilker et all 



12004 IVillanova et al1i2007f) . Still, the sub-populations in uj 
Cen can all be characterised as old and metal-poor. Taking 
uj Cen in this sense as representative for the UCDs, we as- 
sume that their stellar populations are composed of different 
sub-populations, but that these sub-populations are similar 
enough to describe each UCD as a single SSP for the pur- 
pose of this paper. Also note that stellar-encounter-driven 
dynamical evolution is negligible in the UCDs (DHK). A 
disagreement between the SSP models and the observations 
can in this light be interpreted as being caused by assuming 
the wrong IMF. 

The S SP m odels of IBruzual fc Charlotl l|2003l ) and 
IMarastonl (|2005l ) differ by the stellar evolutionary models 
used to calculate the luminosity of the modelled popula- 
tion as well as the total m ass assumed for this population. 
IBruzual fc Charlotl (|2003h assume a somewhat highe r mass- 
loss rate for the stellar populations (|Marastonll2005l . in par- 
ticular her figure 22), while the luminosities they get from 
the stellar models they use are lower. In effect , the esti- 
mates for the M/L ratios by IBruzual fc Charlotl (|2003h are 
similar to the ones by Maraston (2005, her figure 24). How- 
ever, considering the predicti ons for the M/Ly ratios o f old 
populations, the estimates bv lBruzual fc Charlotl J2003I ) arc 
about 20% lower than the ones by IMarastonl (120051 ). Note 
that this cannot be accounted for by the different formu- 
lations for the canonical IMF these authors use since they 



turn out to be nea r ly ide ntical (figure 8 in DHK). In fact, 
IBruzual fc Charlotl |2003h find that the stellar mass of a 
10 Gyr old population is 52% of the initial stellar mass 
for the canonical IMF they use, while it would have been 
54% if th e y had used the same formulation of the IMF as 
IMarastonl (|2005l ) does. The reminder of the difference in 
the M/Ly ratio of an old stellar population must thus be 
the consequence of the different stellar evolutionary mod- 
els used and different assu mptions regarding the remnant 
masses (also see DHK and IMieske et all 2008). As a com- 
promise bet ween the two sets of SSP models, we follow the 
approach by IMieske et al.l ll2008h and take the m ean of the 
predic tions from IBruzual fc Charlotl (J2003) and IMarastonl 
(2005) as the reference for a comparison to the observations 
in UCDs. 

We consider ages of 7 Gyr and 13 Gyr for the UCDs, 
since these values are at the limits of the ages expected for 
them. An age of 7 Gyr would be consistent w ith the inter- 
media te age for the Fornax UCD s suggested in IMieske et al.l 
l|2006t) and IMieske et all l|2008h . Note that assuming even 
younger ages would increase the discrepancy between the ob- 
served M/Lv ratio and the model predictions. Ages higher 
than 13 Gyr are excluded by the estimates fo r the age of the 
universe (13.73to:is Gyr; ISpergel et al.ll2007h . 

The turn-off mass from the main sequence for a popu- 
lation of coeval stars, m to , marks quite well the stellar mass 
above which stars of that population have already evolved 
into stellar remnants. It is ~ 1 Mq for a ~ 10 Gyr stel- 
lar population. Since stellar evolution is slow for old stars, 
m to = 1 Mq is a reasonably good approximation for a 7 Gyr 
old SSP as well as for a 13 Gyr old SSP. 

The contribution of the stellar remnants to the l/-band 
luminosity of the UCDs is small and therefore neglected in 
this paper. The luminosity, L m , of a modelled stellar pop- 
ulation is thus insensitive to the degree of top-heavyness 
of the IMF, since the IMF is only allowed to vary in a 
mass range where the stars have evolved after ~ 10 Gyr. 
The masses of the stars that have not evolved yet are as- 
sumed to be distributed in concordance with the canonical 
IMF. Thus, L m can be determ ined using the SSP m odels 
from IBruzual fc Charlotl (120031 ) and IMarastonl l|2005l ) with 
the canonical IMF. 

We note that by this approach the influence of binary 
systems on stellar evolution is neglected. 



3.1.3 The initial-to- final-mass relation for stars 

In order to find an explicit formulation of M m in eq. (|4]), 
a formulation of the masses of evolved stars as a func- 
tion of their initial masses is needed. This function, called 
the initial-to-final-mass relation, m rom (m), allows to calcu- 
late the total mass of an evolved SSP from its IMF for a 
given age. Using the three-part power-law IMFs from Sec- 
tion 13.1.11 the integral that has to be solved in this calcula- 
tion reads 



TOrem(m)£ p i(m) dm, 



(6) 



where m is the stellar initial mass in M©. The limits of the 
integration are set by the lower and the upper initial mass 
limit for stars. 
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The initial-to-final-mass relation used in this paper is 
specified in the following. 

For stars with initial masses m < mto, m Icm = m is as- 
sumed, i.e. the mass loss of main-sequence stars is neglected. 

Stars with initial masses of m to < m < 8Mq are as- 
sumed to have evolved into white dw arfs (WDs), in con- 
cordan c e with the m a ss lim it given by IKoester fc Reimersl 
(1996). iKalirai et all (|200St ) find, performing a weighted 
least-squares fit of a linear function to data based on ob- 
servations of WDs in star clusters, 

m rcm = (0.109 ± 0.007)^ + (0.394 ± 0.025), (7) 

for a relation between the mass of WDs and the initial mass 
of their progenitors, where m is the stellar initial mass in 
Mq. This relation is adopted in this paper. 

Stars initially more massive than 8 Mq but less massive 
than « 25 Mq are predicted to evolve into neutron stars 
(NSs) with a remarkably narrow mass spread (cf. figures 
12 and 16 in IWooslev et al.ll2002ri. Th is is o bservat ionally 
supported bv lThorsett fc Chakrabartvl (|l999h . who find the 
mass-distribution of pulsars (i.e. observable NSs) in their 
data sample to be consistent with a Gaussian distribution 
with a mean of 1.35 Mq and a width of 0.04 Mq. Thus, in 
this paper 1.35Mq is adopted for the masses of all stellar 
remnants with initial masses between 8 Mq and 25 Mq . 

Stars with initial masses above 25 Mq are generally 
thought to be the progenitors of stellar-mass black holes 
(BHs). However, the theoretical predictions for the masses 
of their remnants are not only strongly dependent on metal- 
licity, but also on the assumptions on how t he evolution of 
such stars proceed s (see figures 12 and 16 in lWooslev et all 
I2002T I. Figure 12 in lWooslev et al ] (|2002r i might suggest that 
the case of the higher remnant masses is the more appropri- 
ate choice for low-metallicity environments such as GCs and 
UCDs. However, the masses of observat ionally confirmed 
BHs lie all in a range that is covered by assuming that the 
remnants of very high-mass st ars only have 10% of the ini- 
tial mass of their progenitors (|Casaresll205"7l y In our paper, 
we thus assume that stars with m > 25 Mq evolve into BHs 
that have either 10% or 50% of the mass of their progenitor 
stars, but the emphasis is on the case with the less massive 
BHs because of the observational support for their existence. 

Note that BHs formed through single-star evolution dif- 
fer from NSs in mass, but not in the processes that precede 
their creation. NSs and BHs are both compact remnants 
that emerge from the core collapse and SN explosion of a 
massive star. 

To summarise, the complete initial-to-final-mass func- 
tion, rrircm used here is 



= < 



ivi ' 

0.109- 
1.35, 
0.1 



+ 0.394, 



M Q < 
mto < 

8€ ^ 



< 25, 



(8) 



25 < wg < m * 



where m to denotes the turn-off mass and m max the upper 
initial mass limit for stars (eqs. [2] and [5}. Inserting eq. ([8]) 
into eq. (J6]l and carrying out the integration on the right 
hand side of eq. (O yields the mass of all stars and stellar 
remnants as a function of only the high- mass IMF-slope, CV3, 
if mto (i.e. age) and m max are specified. The terms resulting 



Table 1. Fit parameters for the metallicity-dependent interpo- 
lation formula for Tv to the data from SS P models with the 
canon ical IMF. BC indicates SS P models fromlBruzual fc CharloiJ 
( 2003) and M SSP models from IMarastorJ l|2005f ). 



Model 


a 


b 


c 


BC, 7 Gyr 


3.29 


0.12 


1.05 


M, 7 Gyr 


3.26 


0.22 


1.24 


BC, 13 Gyr 


3.48 


0.55 


1.71 


M, 13 Gyr 


3.46 


0.79 


1.88 



from this integration for initial stellar masses above m to are 
written down explicitly in Appendix 1X1 



3.1.4 Normalised mass-to-light ratios 

The metallicities estimated for the UCDs usually do not 
coincide with the grid poi nts of the SSP models by 
iBruzual fc CharlotI l|2003u and IMarastonl (|2005h . It is there- 
fore necessary to find interpolation formulae that describe 
the metallicity dependency of the M/Ly ratio in the models 
(which actually is a dependency of the luminosity on metal- 
licity). This can be done by fitting exponential functions of 
the form 



F([Z/H]) 



[Z/Hl+b . 

a' ' ' +c 



Mq 
Lq 



^, 0) 



to the data from IBruzual fc CharlotI l|2003l ) and lMarasto"nl 
(2005), where [Z/R] is the metallicity. The best-fitting pa- 
rameters a, b and c found in a least-squares fit to the mod- 
els used in this paper are listed in Table [1] The excel- 
lent agreement of this typ e of function to the m odels from 
IBruzual fc CharlotI (|2003l ) and IMarastonl (|2005l ) is demon- 
strated in fig. 9 in DHK. 

The reference relation that is taken to describe the 
metallicity dependency of the M/Ly ratio for a SSP with 
a certain age and with the canonical IMF is the mean of 
the c orresponding relations der ived from the SSP m odels 
from | Bruzual fc CharlotI i|2003D and IMarastonl i|2005h (cf. 
iMieske et alj|2008h . The ratio between the observed M/L v 
ratio for a UCD and the result from the reference relation at 
the appropriate metallicity is a measure for the discrepancy 
between the observed value and the theoretical prediction. 
It is convenient for the purpose here to multiply these val- 
ues by the prediction of the reference relation for the M/Lv 
ratio at Solar metallicity. These quantities will be referred 
to as normalised M/Ly ratios, Ty in , 



Fbc\z + Fm\z 



x (F BC \z e +F m \z q ), (10) 



where a subscript BC indicates t hat the parameters a , b and 
c correspond to a SSP model from lBruzual fc CharlotI (I2003T ) 
and the subscript M indicates tha t the parameters a, b and 
c correspond to a SSP model from Maraston (2005) (for the 
same age). 

Using these values for Ty,n, eq. Q can be rewritten as 



T v ,n = 0. 



(11) 



Lm\z Q is thereby no longer a metallicity-dependent variable, 
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but is fixed to the value the reference relation predicts for So- 
lar metallicity and thereby only depends on the age assumed 
in the model and the amplitude of the factors fc in the IMF. 
The metallicity dependency is shifted into the transforma- 
tion from the observed M/Lv ratio of the UCD to Ty, n . 
The Tv,n values are noted in Tab. and shown in Fig. [T] 
Their uncertainties have been propagated from the errors of 
the observed dynamical M/Lv ratios and the errors of the 
metallicity estimates. 

The numerical value of L m is calculated from the sec- 
ondary condition that the prediction for Tv.n from the SSP 
models should correspond to a stellar population with the 
canonical IMF with m max = 100 Mq and a full population 
of remnants (i.e. a stellar population as in the SSP models). 
For this, m rcm as given in eq. © is used, adopting the case 
that the black hole masses, tubh, are 10% the stellar initial 
masses (ttibh = 0.1m). 

There is evide nce that GCs u s ually have q-enrichments , 
[a/Fe], of 0.3 dex (|Carnevlll996h . lEvstigneeva et all \200H \ 
find that the same [a/Fe] is also typical for the UCDs in the 
Virgo cluster they examine. On the other hand. lMieske et all 
|2007l ) find that a number of UCD can didates is consistent 
with having Solar [a/Fe], which is why iMieske et all HH) 
adopt Solar [a/Fe] for all stellar systems in their study. How- 
ever, assuming a super Solar [a/Fe] is the more careful choice 
in the context of the present paper, since it attributes more 
of a possibly enhanced M/Lv ratio in UCDs to metallicity 
effects. As in DHK, we therefore adopt [a/Fe] = 0.3 dex for 
all GCs and UCDs and estimate their metallicities, [Z/B], 
from their iron abundances, [Fe/H]. This is done using the 
relation 

[Z/B] = [Fe/H] + 0.94 [a/Fe] (12) 

taken from I Thomas et al.l (|2003r ). Consequently, the [Z/H] 
used for calculating the Tv.n are 0.28 dex higher than the 
[Fe/H] and the Ty, n in this paper ar e thereby slightly lower 
than the ones in lMieske et alj (|2008l '). 

The assumed age turns out to be almost irrelevant for 
the Tv,n calculated for the individual stellar systems. How- 
ever, the assumed age does have a strong impact on the Tv,n 
predicted by the SSP- models (see also DHK). 



3.1.5 The fate of the processed material and the stellar 
remnants in the UCDs 

In order to have an influence on the dynamics of a stellar 
system, the stellar remnants that form in it have to remain 
bound to it. This can be assumed to be the case for the WDs 
in the UCDs, since WDs inherit the peculiar velocities of 
their progenitor stars and two-body encounter driven mass 
loss is negligible for the UCDs (see Section [37L2]) . 

Unlike the case with WDs, stellar evolution has a direct 
impact on the velocity distribution of NSs. It is well estab- 
lished that many pulsars move with high peculiar velocities, 
which they must have obtai ned somehow in their formation 
out o f their progenitor stars jWooslevll 19871 ; iLvne fc Lorimerl 
ll994T ). lLvne fc Lorimerl (1 19941 ) give the mean pulsar birth ve- 
locity as 450 ± 90 km s~ . Since the processes that lead to 
the formation of BHs through single-star evolution are the 
same as the ones that precede the formation of NSs, the BHs 
should also receive kicks. 



The UCDs have velocity dispersions of < 50km s _1 , 
which suggests escape velocities of the order of < 100km s _1 . 
Thus, the peculiar velocities of most NSs and BHs should be 
high enough to leave the UCDs. On the other hand, NSs are 
known to populate GCs, which suggests that also the UCDs 
are able to retain some fraction of these objects. 

Most of the matter processed in intermediate- mass and 
high-mass stars is reinserted as gas and dust into the inter- 
stellar medium during stellar evolution. Its fate is therefore 
crucial for the developement and consequently the M/Lv 
ratio of a stellar system. 

There are in general three possibilities for what can hap- 
pen to this material. If it remains inside the cluster, it can 
(at least in principle) simply accumulate (and thereby emit 
almost no radiation in the V-band) or it can be used up 
in the formation of subsequent stellar populations. Alter- 
natively, the gas can be driven out of the UCDs, e.g. by 
type I SNe or by the ram pressure caused by the movement 
of the UCD through the intergalactic medium. 

Gas and dust originating from intermediate-mass stars 
has a good chance to stay inside the UCDs, since these 
stars form in their final stage planetary nebula e that expand 
with moderate velocities (as 20kms _1 ; see e.g. lGesicki et all 
2003). These velocities are too low for the matter to leave 
a star cluster with a deep potential well immediately. This 
makes massive AGB stars attractive pr ogenitors for a sec- 
ond g enera tion of stars, as propo sed in iD'Antona fc Caloil 
J200 



generation ot stars, as propos 
2) and lD'Antona et all \200H \. 



However, massive stars evolve into SN and thereby re- 
lease as 10 50 erg per Sol ar unit of initial ma ss of the pro- 
genitor star (cf. fig. 1 in lNomoto et al.l 120061 ) . This clearly 
exceeds the binding energy of a star to a UCD. Mate- 
rial originating from these stars will therefore easily escape 
from the stellar system, unless the kinetic energy of the gas 
from the SN explosion is dissipated (e.g. by the interaction 
with primordial gas or the collision of expanding gas en- 
velopes from different SNe with one another). The gas den- 
sity and the holding time inside the UCDs might become 
lo ng enough for the gas to co ol and to collapse, as discussed 
in lTenorio-Tagle et al.l < |2007h . 

Note however that neither self-enrichment by massive 
AGB stars nor self-enrichment with SN ejecta can explain 
the multiple stellar populations in uj Cen, since both sce- 
narios act on a time scale of < 200 Myr, whereas the age 
difference betw een the different stellar populations in i o Cen 
is a few Gyr 



IVillanova et al.l 



Hilk er fc Richtl er 2000; [Hilker et al.l [2004; 



2007) 



The essence of this is that the current knowledge on 
the evolution of the UCDs does not allow solid conclusions 
on the composition of the UCDs. We therefore consider six 
different compositions of the UCDs for which we estimate 
the high-mass IMF-slope: 

(i) Out of all material from burnt out stars, only WDs are 
retained by the UCDs. This can be taken as the lower limit 
for the amount of matter that stays inside the UCDs since 
the UCDs are nearly unaffected by dynamical evolution (cf. 
Section [3X2]) . 

(ii) 20% of the compact remnants from stars initially 
more massive than 8M© are retained by the UCDs. The 
remnant masses of stars with m > 25 Mq, ttibh, are as- 
sumed to be 10% of the initial mass of their progenitors, 
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Figure 1. Normalised M/Ly ratios, Ty n , of the stellar systems collected in table 5 in lMieske fc Kroupal l|2008h . provided their dynamical 
mass is estimated to be 2 X 10 6 Mq or more. Thus, the figure only shows objects that are UCDs according to the definition used in this 
paper. The ages assumed for them are either 13 Gyr (left panel) or 7 Gyr (right panel). The dashed horizontal lines indicate Yy a for a 
SSP that formed with the canonical IMF and is of the age that is assumed for the UCDs in the according panels. The dotted horizontal 
lines correspond to the mean of the Tv,n of all UCDs in the sample and the shaded areas indicate the uncertainty given to this value. 
These numbers are used to estimate the high mass-slope of the UCDs (see Section l4.2.1H . 



771bh = 0.1m. The NS and BH retention rate of 20% is an 
arbitrarily chosen value, but this scenario might still be close 
to a realistic one. On the one hand it allows for some NSs 
in the UCDs as observed for GCs, but on the other hand it 
also takes into account that the observed velocity dispersion 
of pulsars is high by estimating the fraction of retained NSs 
to be low. 

(iii) As scenario (ii), but with ttibh = 0.5m. 

(iv) All stellar remnants are retained by the UCDs and 
ibbh = 0.1m. Such a population, where all stars and stellar 
remnants, but not the ejecta from stars are considered, is 
assumed in the SSP models. 

(v) As scenario (iv), but with itibh = 0.5m. 

(vi) The UCDs were gas-free after star-formation ceased 
in them, but all material that was processed in burnt-out 
stars is retained by the UCDs and star formation with the 
gaseous component of this matter is somehow inhibited. Of 
all the models considered here, this is the one where stars 
contribute the least to the total mass of the UCD (con- 
sisting of stars, remnants and possibly gas). Note however 
that t he hydrodynamic calculations by iTenorio-Tagle et al.l 
(|2007l ) suggest that such a scenario is unlikely because the 
gas accumulating in the UCDs due to stellar evolution will 
more likely either leave the UCDs or collapse into new stars. 

Interstellar gas and all remnants in the UCDs are con- 
sidered not to contribute to the light of the UCDs. In other 



words, the very high M/Lv ratios of these components of 
the UCDs are taken to be infinity. 



4 RESULTS 

The value for the high-mass IMF slope, 0:3, implied by a 
given normalised M/Lv ratio, Tv,n, for the assumptions on 
the stellar populations in the UCDs specified in Section [3] 
can be calculated from eq. {TTJ using the Newton-Raphson 
root-finding method. 

There is a lower limit for the Ty,n that leads to a so- 
lution for eq. |TT}, because the lowest Tv,n that can be 
realised within the model is the one for a stellar population 
whose IMF is cut off at 1 Mq (i.e. £13 = 00). If the age of 
the UCDs is assumed to be 13 Gyr, the individual Ty,n of 
a number of UCDs is actually below that limit. 

Close to that lower limit, «3 increases rapidly with de- 
creasing Tv,n, implying a steep high-mass IMF (Fig. [2}. In 
this range, the solutions to eq. (|lip become degenerate for 
the different assumptions on how much mass is retained by 
the UCDs. That is because a steep high-mass IMF means 
few high-mass stars and it is therefore not decisive for the 
Tv,n of an old stellar system how much matter from those 
stars is retained. 

The CV3 of the canonical IMF is in the regime where 
the relation between Tv> and 03 is already close to be- 
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Figure 2. The high-mass slope, 0:3, against the normalised mass-to-light ratio, Tp„, as implied by solving eq. for a three-part 

power-law IMF (see eq. [5]and the equations in the appendix). The curves are for a 13 Gyr old SSP (upper panel) and for a 7 Gyr old SSP 
(lower panel). The different styles of the curves correspond to different assumptions on how much processed matter (with extremely high 
M/Ly ratio) besides WDs is retained by the UCDs (from the bottom to the top curves in each panel): no remnants of massive stars; 
20% of the remnants of massive stars; all remnants of massive stars; all material processed by burnt-out stars. Two curves of the same 
style indicate different assumptions for the upper mass limit of the IMF for the same assumption on the matter retained in the UCDs: 
100 Mq (lower curve) and 150 Mq (upper curve). The dashed horizontal lines indicate in each panel the canonical high-mass IMF index, 
03 = 2.3. Its intersections with the curves show the Ty,n which the canonical IMF would imply for a particular remnant population. 
The Ty n of the individual UCDs are shown as crosses at the bottom of each panel. The dashed vertical lines indicate Ty n for a SSP 
that formed with the canonical IMF and is of the age that is assumed for the UCDs in the according panels. The dotted vertical lines 
correspond to the mean of the Ty }11 of all UCDs in the sample and the shaded areas indicate the uncertainty given to this value (see 
Section 14.2. It . The intersections of a vertical line with the curves show that 03 corresponding to a particular Ty „ for the different 
assumptions on the retained remnant population. In this Figure, the remnants of stars with m > 25 Mq are assumed to have masses of 
10% of the initial mass of their progenitors whereever this is relevant. 



ing degenerate for different assumptions on the NS and 
BH retention rate. Only assuming all matter processed in 
burnt-out stars remains inside the UCDs without forming 
new stars would lead to a distinctively higher TV, n for 
03 = 2.3. In other words, the predictions of the models from 
iBruzual fc Chariot! (|2003l ) and lMarastonl |2005h for the TV ,n 
for a SSP with the canonical IMF (shown as the horizontal 
dashed line in Fig.[T]and as the vertical dashed line in Fig. [3]) 
depend, except for extreme cases, only weakly on the fate 
assumed for the material processed in massive stars. 

Comparing solutions of eq. (|11[) for the same YV,n and 



remnant retention rate, but for upper mass limits of 100 Mq 
and 150 Mq, reveals that the remnants of very massive stars 
do not play a decisive role for the «3 that are obtained, as 
illustrated in Figs. [2] and [4] This finding may be surpris- 
ing, since the total mass of the remnants of high-mass stars 
is a function of the exponent Q3 (see Appendix \K§ . This 
mass must therefore increase dramatically with increasing 
Q3 above some critical value for 03. 

The results of solving eq. (|lip if a remnant retention 
rate of 20% and ?tibh = 0.1m is assumed are noted in Ta- 
ble [2] However, for many individual UCDs solutions do not 
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exist if a high age is assumed for them (i.e. their Tv,n 
is clearly below the prediction from the SSP models for 
a canonical IMF), and the uncertainties are large in any 
case. On the other han d, application of the Pearson tes t for 
the goodness of fit (cf. iBhattacharvva fc Johnspnlll977l and 
DHK) on the 46 UCDs from iMieske et all IMPOST ) shows that 
the actual distribution of the Tv,n of the UCDs in the sam- 
ple is highly unlikely if their individual Tv,n scatters equally 
to both sides of the prediction for the Ty,, of a SSP with 
the canonical IMF (less than 1% if the age of the UCDs 
is assumed to be at its maximum, 13 Gyr, and much less 
than 0.5% if the age of the UCDs is assumed to be 7 Gyr). 
The properties of the sample of UCDs therefore imply an 
IMF that deviates from the canonical IMF (provided that 
they do not contain non-baryonic DM), such as a three-part 
power-law IMF with Q3 < 2.3. 

The emphasis in this paper is therefore on constraining 
likely values for the high-mass IMF slopes of the UCDs from 
the properties of the whole sample of UCDs and different 
subsamples thereof. It is decisive for this to know whether 
the Yy,n of the UCDs are correlated with their mass, M, 
and to quantify this correlation if there is one (Section 14. 
If such a dependency is found, the dependency of Q3 on Ty. n 
can be translated into a dependency of 03 on M. 



4.1 Does «3 depend on mass? 

At present, it is unclear whether UCD s are the most mas - 
sive GCs (e.g. IMieske et alj|2002l , liooi iForbes et alj|2008h 
or whether UCDs and GCs are different p opulations (e.g. 
iDrinkwater et al.l I200I iGoerdt et al] 120081 ). However, the 
answer to this question has implications on how a depen- 
dency of Tv,n on M, Tv,n{M), has to be formulated for 
GCs and UCDs. An appropriate formulation of Yy-, n (M) as 
a representation for the typical Yy*, n of objects with a given 
mass would be a single, continuous function in the first case, 
but different functions for GCs and UCDs in the second case. 

The MWGCs, which make up most of th e GCs in the 
sample used here (tab. 5 in IMieske et al.ll200Sf ) , show no ev- 
idenc e for a bulk-dependency of Ty with M (|McLaughlinl 
2000). Therefore, the mean Yy-, n of the GCs in the data sam- 
ple, Tgc, is adopted for Yv, n (M) in this mass range. Thus, 
for M < 2 x 10 6 M Q , Tv>(M) = 2.43 ± 0.16 M h~] v if the 
assumed age is 7 Gyr and T v ,n(M) = 2.61 ± 0.18 M L" 1 ^ 
if the assumed age is 13 Gyr (uncertainties are one-sigma 
values) . 

The uncertainties of the data for the UCDs leaves many 
options for an appropriate formulation of Tv,n(M) for them. 
We choose 



Tv,n(M) 



A 



lop 



+ B 



M© 
L 
(13) 

for M > 2 x 10 6 M Q , where M is in Solar units and A and B 
are parameters which are either fixed by a secondary condi- 
tion or determined by a least-squares fit. Note that weight- 
ing the uncertainties when fitting is not advisable in this 
case, as it would cause an unwanted bias. This becomes evi- 
dent by considering two stellar systems with the same mass 
and uncertainty of the mass, but different luminosities. The 
uncertainty of a luminosity measurement is negligible com- 
pared to the uncertainty of a mass estimate. The uncertainty 



of the M/L ratio is thus higher for the stellar system with 
the higher M/L ratio, even if the parameter that induces 
this uncertainty is the same for both systems. The param- 
eters A and B are therefore determined with equal weight 
to every measurement and the uncertainties of A and B are 
estimated only from the scatter of the data. 

In order to constrain Tvn(M) for the UCDs in the case 
that UCDs and GCs are two distinct populations, A and B 
in eq. (|13|l are left as free parameters for the fit. The best- 
fitting parameters are A = 1.84±0.89 and B = 3.71 ±0.70 if 
the UCDs are assumed to be 7 Gyr old and A = 1.87 ± 0.81 
and B = 3.47 ± 0.64 if the UCDs are assumed to be 13 Gyr 
old. This may hint at a systematic increase of the Tv,n of 
the UCDs with M, but the significance of this result (« 2a) 
is not high enough to allow definite concl usionj]]. This find- 
ing is consistent with IMieske et all (|2008l ). who performed a 
similar test but only for the UCDs in Fornax. 

If UCDs are the most massive GCs, Ty jn (M) is ex- 
pected to be continuous at M — 2 x 10 6 M , which in this 
paper is taken to be the mass that separates GCs and UCDs 
(see Section[2j|. For this case, Tv,n{M) of the UCDs is there- 
fore estimated by setting B in eq. (|13p to the numerical value 
of Tgc in Solar units and leaving only A as a free parameter 
to be determined in the fit. The result is A = 3.00 ± 0.42 if 
the UCDs are assumed to be 7 Gyr old and A = 3.04 ± 0.46 
if the UCDs are assumed to be 13 Gyr old. In this case, the 
increase of the Tv, n of the UCDs with their mass is highly 
significant. 

Tv,n(M) is plotted in Fig. [3] together with the data for 
the GCs and the UCDs. 



4.2 Constraining 03 from the whole sample of 
UCDs 

4-2.1 UCDs and GCs as independent populations 



It was shown in Section 14.11 that there is no hard evidence 
for a correlation of the Ty jn of the UCDs with their mass if 
UCDs and GCs are separate populations. For this case, it is 
therefore a useful and good assumption that all UCDs have 
the same Tv,n and that deviations from it are due to statis- 
tical scatter. Thus, we estimate the Ty-, n of the UCDs and 
the uncertainty of this value by performing a least-squares 
fit of eq. {I3J with A = to the T v ,n of the UCDs. The 
best-fitting parameter B then equals to the numerical value 
of the mean of the Yy,n of the individual UCDs, Tucd, and 
is (4.75 ± 0.34) M L~* v if the age of the UCDs is assumed 
to be 7 Gyr and (4.97±0.37) M L~^ v if the age of the UCDs 
is assumed to be 13 Gyr. Tucd is shown as the horizontal 
dotted line in Fig. [1] and as the vertical dotted line in Fig. [2] 
The uncertainties of Tucd are indicated as shaded areas in 
these figures. 

A given value for Tucd is taken to depend only on 03. 
This implies that all UCDs have formed with the same (top- 
heavy) IMF, which can be considered as characteristic for 
very dense star- forming regions. Note however that even if 
this assumption is consistent with the available data, it is 



1 We mention that weighting the UCDs by the uncertainties leads 
qualitatively to the same results, although the best-fitting values 
for A and B und their uncertainties are slightly lower. 
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Table 2. Normalised M/Ly ratios, Ty^, of the UCDs for assumed ages of 7Gyr and 13 Gyr, and the high-mass slopes, 03, these TV, n 
suggest if 20 % of the remnants of massive stars are retained by the UCDs, BHs have 10% of the initial mass of their progenitor stars and 
the u pper mass limit of the IMF is m max = 100 Mq. The contents of the columns are the following: Column 1: The object identification 
(as in lMieske et al,l|2008l . table 5), Column 2: The projected half-light radius of the UCD, Column 3: The esti mate for the iron-abun dance, 
Column 4: The mass of the UCD, Column 5: Its Ty n based on the models bv lBruzual Sz C hariot ( 2 003! ) and iMarastonl {2005) for a 
7 Gyr old SSP with the canonical IMF, Column 6: The estimate for 03 based on the value for Ty.a in Column 5, Columns 7 and 8: As 
Columns 5 and 6 respe ctively, but fo r an as su med age of the UCD o f 13 Gyr. The superscript nu mbers in Colu mn 3 indicate the origin of 
the [Fo/H] estimate: b lMieske et al.l (120081 1. 2: lHasegan et al.l <2005ft . 3: lMevIan et al.l teOOlT ), 4: lHarrisl fl996). A superscript * indicates 
that [Fe/H] was not obtained from colour indices, but from line indices or the properties of the resolved stellar population of the stellar 
object (private communication with S. Mieske). Ty.n is estimated using [Z/H], which is £s 0.3 dex higher than the corresponding [Fe/H] 
due to the assumed o-enhancement of the stellar systems (see Section |3. 1.411 . Dots in Columns 5 and 8 indicate where no solution for 
eq. I|ll[l is found under the given assumptions. The canonical IMF would have 012 = 0:3 = 2.3 (Salpeter-Massey index, eq.[2j. 
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Figure 3. Normalised M/Ly ratios, Ty in , of all objects listed in table 5 in lMieske et al.l l|2008h . The assumed ages for them are either 
13 Gyr (left panel) or 7 Gyr (right panel). Contrary to Fig. [T] this figure also shows stellar systems with dynamical masses less than 
2 X 10 6 Mq, i.e. stellar systems that are not UCDs but GCs according to the definition used in this paper. Objects considered as GCs are 
scpcratcd from objects considered as UCDs by the thin, dashed vertical line in each panel. The solid line indicates Yy n (M), a function 
that describes the systematic increase of the average Tv, n with mass for the case that GCs and UCDs are a single population. If GCs 
and UCDs are separate populations, Yy n (M) of the UCDs is represented by the dotted line that starts at 2 X 10 6 Mq. Note that in 
this case the uncertainty to the slope is very high and is therefore not significant. The dashed horizontal line indicates the prediction for 
T\/. n of an SSP with the canonical IMF. 



a simplification because the TV, n of the UCDs are expected 
to scatter due to age differences. Furthermore, the sugges- 
tion that the IMF is top-heavy in UCDs in comparison to 
less massive stellar systems is based on the notion that the 
process of star formation (and thus the IMF) depends on 
the physical conditions under which it takes place. This im- 
plies that there can be as many IMFs as physical conditions 
under which star formation takes place. 

Tucd can be translated into different expected values 
for ay, depending on the assumed remnant population in the 
UCDs. The upper limit to the expected as can be obtained 
from the lower bound of the uncertainty of Tucd, and the 
lower limit to the expected 03 can be obtained from the 
upper bound of the uncertainty of Tucd. The values for «3 
are listed in Tab. [3] 

iMieske et aLI (|2006l ) found that the UCDs in Fornax 
and the ones in Virgo are (despite their similarity) distinct 
in their properties, which could indicate a different origin or 
age for the two groups. Among these distinctive properties 
is the average Yy <n of the UCDs, which is clearly higher for 
the ones in the Virgo Cluster. It is therefore worthwhile to 
relax the assumption of a common Tv.n for all UCDs and 
assume a common Tv,n only for the UCDs that are bound to 
the same larger structure (i.e. the Fornax Cluster, the Virgo 
Cluster or Centaurus A), although the smaller size of the 



subsamples decreases their statistical significance. Estimates 
for Q3 for the UCDs in the different subsamples can then 
be obtained in the same way as for the whole sample. The 
results are given in Tab. [3] 

Note that the values obtained for 03 for the UCDs in 
Fornax have large uncertainties if they are assumed to be 
13 Gyr old. This is because in that case Tucd for them is 
quite close to the Ty in where ctz asympotically approaches 
infinity. Consequently Q3 is obtained from an interval in 
Tv, n where small variations of Ty, n imply large changes 
in 03. This applies in particular to the upper limit to a.3. 
Thus, the numbers in question are only of use for giving 
lower bounds for the high-mass IMF slope, which are ob- 
tained at a T^, n where the dependency of «3 on Tv, n is 
more moderate. 

Also note the strong impact of S999 with its extreme 
T^,n (upper-most data-point in Figs. [T] and |3jon the Q3 
derived for the UCDs in the Virgo Cluster. The relevance of 
this particular cluster is evidently much smaller if the whole 
sample of UCDs is considered. The main results presented 
in this paper are therefore either not or only mildly affected 
by this stellar system. In particular, S999 plays no role for 
deciding whether there is a significant correlation between 
the Tv, n of the UCDs and their mass. Such an outlier may 
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Table 3. Estimates for the most likely values of az for different assumptions concerning age and remnant population of the UCDs, as 
detailed in Section 14,2.11 The first column specifics the supposed remnant population of the UCDs (m max = 100 Mq). For the stars more 
massive than 25 Mq, the cases of them forming BHs with 10% their initial mass, itibh = 0.1m, or 50% their initial mass. rrtBH = 0.5m, 
are considered whereever this makes a difference. The second column displays the 03 corresponding to the mean of the Ty n of the UCDs 
and the uncertainties to 03 calculated from the uncertainties to the mean Ty n . This can be taken as a convenient number to quantify 
the high-mass IMF slope of the UCDs as a class of objects. The numbers in Columns 3 to 5 have the same meaning as the numbers in 
Column 2, but for different subsamples of UCDs. The subsamples are chosen by the larger structures the UCDs are bound to, namely 
the Fornax Cluster (Column 3), the Virgo Cluster with S999 (Column 4) and without S999 (Column 5) and Centaurus A (Column 6). 





All Fornax Virgo 


Virgo 


Centaurus A 




(with S999) 


(without S999) 




Model 


03 03 S3 


S3 


«3 


assumed age of 13 Gyr 



no SN remnants 
20% of the SN remnants, meH 
20% of the SN remnants, ttibh 
all SN remnants, meH = 0.1m 
all SN remnants, rngj-i = 0.5m 
all processed material 



: 0.1m 
: 0.5m 



1.35 
1.57 
1.78 
1.85 
2.11 
2.93 



0.23 
0.17 
+0.17 
0.12 
+0.13 
0.09 
+0.14 
0.10 
+0.11 
0.08 
+0.26 
0.18 



2.10 
2.17 
2.26 
2.33 
2.50 
3.90 



0.90 
0.42 
+0.84 
0.35 
+0.76 
0.29 
+0.73 
0.29 
+0.01 
0.23 
1.55 
0.59 



0.81 
1.22 
1.50 
1.56 
1.86 
2.43 



+0.37 
0.23 
+0.23 
0.13 
+0.19 
0.11 
+0.20 
0.12 
+0.18 
0.10 
+0.32 
0.18 



1.11 
1.41 

1.65 
1.72 
2.00 
2.69 



+0.25 
0.18 
+0.17 
0.12 
+0.14 
0.09 
+0.14 
0.10 
+0.12 
0.08 
+0.25 
0.16 



1.49 
1 

1.86 
1.94 
2.18 

3.08 



0.29 
0.21 
+0.23 
-0.15 
+0.18 
0.12 
+0.18 
0.12 
+ 0.15 
0.10 
+0.36 
0.22 



assumed age of 7 Gyr 



no SN remnants 
20% of the SN remnants, rajn 
20% of the SN remnants, rajn 
all SN remnants, mgH = 0.1m 
all SN remnants, meH = 0.5m 
all processed material 



: 0.1m 
: 0.5m 



4G+°'09 
u, ^ a -0.08 

1 04+ ' 05 

1-0.04 
-0.04 
+0.05 
0.04 
+0.04 
0.04 
+0.07 
0.05 



1.34 
1.40 
1.72 
2.19 



0.73 
1.17 
1.46 
1.52 
1.82 
2.36 



+0.15 
0.12 
+0.09 
0.07 
+0.07 
0.06 
+0.08 
0.06 
+0.07 
0.05 
+0.12 
0.09 



0.19 
0.88 
1.21 
1.25 
1.59 
2.00 



+0.20 

0.15 
+0.10 
0.08 
+0.09 
0.07 
+0.10 
0.07 
+0.09 
0.06 
+0.13 
0.09 



0.36 
0.97 
1.28 
1.33 
1.66 
2.10 



+0.11 
0.10 

+0.07 
0.05 
+0.05 
0.04 
+0.06 
0.05 
+0.05 
0.04 
+0.08 
0.06 



0.56 
1.09 
1.38 
1.43 
1.75 
2.25 



+0.11 

0.09 
+0.06 
0.05 
+0.04 
0.04 
+0.05 
0.05 
+0.05 
0.04 
+ 0.08 
0.07 



be du e to recently induced tidal effects (|Fellhauer fc Kroupal 

l2ooeh . 



4-2.2 GCs and UCDs as a single population 

Contrary to the case that UCDs and GCs constitute different 
populations, the slope A is highly significant if UCDs and 
GCs are a single population. A relation between 03 and 
mass, 03 (M), can be established by solving eq. (|lip for the 
different TV,n(M) corresponding to different masses. The 
results for this are plotted in Fig. [4] The fact that in the mass 
range of GCs, Tv, n (M) almost coincides with the model 
predictions for a SSP with the canonical IMF for an assumed 
age of 7 Gyr is purely coincidental. Independent estimates 
on the ages of MWGCs (which are the bulk of the GCs in 
the sample used here) are closer to 13 Gyr than to 7 Gy r 
for most of them (|VandenBerdl2000l ; ISalaris fc Weissll2002t) . 
However, if an age of 13 Gyr is assumed for the GCs and 
UCDs, a 3 (M) is not defined for M < 3 x 10 6 M Q because 
in this mass range TV, a (M) is below the minimum value for 
which eq. (|11|) is solvable. 

This finding implies that the present-day stellar mass 
function of the corresponding stellar systems has to be 
poorer in very low-mass stars than the canonical IMF, since 
their Yy jn cannot be realised in any case if the mass func- 
tion of their main sequence stars equals the canonical IMF. 
There are different processes which tend to drive very low- 
mass stars out of a star-cluster. One of them is dynami- 
cal e volution (cf. lKruiisse"nl [2008 and Kruiisse n fc Lamersl 
2008). It acts faster the less massive a star cluster is. How- 



ever, the exp ected effect on the M / Ly ra tio is only small 
acc ording to | Baumgardt fc M akino ( 20"03l) . their figure 14 



and lBorch et al.l (|2007f ). More relevant would be gas expul- 
sion if the GCs were initially mass-segregated ( [Marks et al.l 
2008). In other words, if GCs have formed with the canon- 
ical IMF, their stellar mass f unctions must hav e changed 
with time (see also DHK and lMieske et all I200SI ). The de- 
pendenc y of the stellar mass fu nction on the cluster concen- 
tration |De Marchi et all 120071 ) indeed suggests this to be 
the case. Consequently, the assumption of a stellar popula- 
tion only altered by stellar evolution would only be valid for 
stellar systems with M > 3 X 10 6 M . 

For a parametrisation of 03 as a function of the mass 
of a stellar system we suggest 



as(M) 



logx 



0.85(1Q- 6 M/M Q ) 

(io~ 6 m/Mq) - : 



+ 0.42, (14) 



for M > 2 x 10 6 M , and oi'(M) = 2.3 for M < 2 X 10 B M 0) 
where M is measured in M . of (M) thus returns the canon- 
ical IMF for GCs, which is motivated w ith the invariance of 
the IMF in resolve d stellar populations (|Kroupall200ll , 120081 ; 
Marks et al.l I2OO8I ). In the range of massive UCDs, which 
are the least vulnerable to dynamical evolution, ct3~(M) is 
chosen to be roughly the mean of 03 (M) for assumed ages 
of 7 Gyr and 13 Gyr at a NS and BH retention rate of 20%. 
Note that the change of as(M) in this mass range is only 
moderate for the two extreme assumptions on the age of 
the UCDs. In the intermediate mass range from 2 X 10 6 M© 
tO « 1O 7 M , 03~(M) is an (in principle arbitrary) interpola- 
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Figure 4. The high-mass slope, 03, in dependency of mass M, as implied by solving eq. l|ll|l for the values of Yy- !n (M) at the according 
M and for a three-part power-law IMF (see eq. [5] and the equations in the Appendix). The curves are for a 13 Gyr old SSP (upper panel) 
and for a 7 Gyr old SSP (lower panel). The different styles of the curves correspond to different assumptions on how much processed 
matter (with very high M/Ly ratio) besides WDs is retained by the UCDs (from the bottom to the top curves in each panel): no 
remnants of massive stars; 20% of the remnants of massive stars; all remnants of massive stars; all material processed by burnt-out stars. 
Two lines of the same style indicate different assumptions for the upper mass limit of the IMF for the same assumption on the matter 
retained in the UCDs: 100 Mq (lower line) and 150 Mq (upper line). The dashed horizontal lines indicate in each panel the canonical 
high-mass IMF index, 0:3 = 2.3. In this Figure, the remnants of stars with m > 25 Mq are assumed to have masses of 10% of the initial 
mass of their progenitors whereever this is relevant. 



tion from the low-mass regime to the very high-mass regime. 
Q3~(Af) is plotted in Fig. [5] 



5 DISCUSSION 

5.1 Stability of the UCDs 

iBaumgardt et all (|2008l ) show that for a star cluster with the 
canonical IMF the combined energy input from all SNe ex- 
ceeds the binding energy of star clusters with initial masses 
up to « 10 7 Mq (cf. their fig. 3 j3- This implies that star 



clusters loose not only most of the matter bound in mas- 
sive stars, but also their primordial gas in less than 40Myr 
(« 40Myr is the time it takes un til all massive stars in a 
SSP have evolved, cf. the grids bv lSchaller et~aIlll992D . Al- 
though many UCDs certainly had initial masses higher than 
10 7 Mq , they also had many more massive stars that evolved 
into SNe if they formed with IMFs as top-heavy as sug- 
gested in this paper. In this case, they would loose an even 
larger fraction of their initial mass during their early evolu- 
tion than less massive stellar systems, because of the large 
mass-faction bound in massive stars. For instance, 23.0% of 



2 Star clusters less massive than 10 7 M© can survive this energy 



input because the energy from the SNe is not distributed uni- 
formly on all matter in the cluster. 
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Figure 5. The parametrisation of the high mass IMF slope, 03, 
as a function of the mass of a stellar system, M, as given by 
eq. 1141) . It is indicated by the solid line. Also shown are a^{M) 
as found from eq. {TTJ using eq. I|13|l assuming a NS and BH 
retention rate of 20% and BHs having 10% of the initial mass of 
their progenitor stars. The upper dashed line corresponds to an 
estimated age of 13 Gyr and the lower dashed line to an estimated 
age of 7 Gyr for the UCDs. 



the total initial stellar mass of a star cluster is in stars more 
massive than 8M© if Q3 = 2.3. This value rises to 73.0% 
for 03 = 1.57 (which is the high-mass IMF slope suggested 
in Tab.[3]for 13 Gyr old UCDs that retain 20% of their NSs 
and BHs with ttibh = 0.1m) and to 93.3% for 0:3 = 1.04 
(which is the corresponding value for 7 Gyr old UCDs). 

Observations of star-forming regions in the Milky Way 
show that only a frac tion of available gas is actually con- 
verted into stars (e.g. lLada fc Ladal [20031 ) . Assuming that 
this left-over gas is swept out of young star clusters by the 
radiation and evolution of massive stars, the total mass-loss 
of the stellar system until the end of massive-star evolution 
can be written as 



Mfi^ = Mtait[l-SFE(l- a:) 



(15) 



where Mj n it and Af nna i are the initial mass of the stellar 
system and the final mass of the stellar system respectively 
(stars and gas), SFE is star formation efficiency of the stellar 
system and x the mass- fraction of stars with m > 8 Mq . It is 
thereby assumed in eq. (|15l) that the total mass of the rem- 
nants remaining in the stellar system is negligible compared 
to the total initial mass of all stars with m > 8M©. This 
approximation is well fulfilled for all cases where the mass of 
the remnants of massive stars is rather small or only a few 
of them remain inside the stellar system. (These cases corre- 
spond to our models i, ii, iii and iv for the composition of the 
UCDs; see Section ^. 1.51 That is why the IMFs estimated for 
the UCDs turn out to be so fiat if one of these models is as- 
sumed for their composition, even if the difference between 
the predicted M/L ratio and the observed M/L ratio is not 
very large.) Assuming that the star formation efficiency in 
a UCD is 0.4, the mass-loss within the first 40Myr is, ac- 
cording to eq. 1151 69.2% of the total initial mass (stellar and 
gas) if the high-mass IMF slope of the UCD was 03 = 2.3 



(canonical IMF), but 89.2% for a 3 = 1.57 and 97.3% for 
Q3 = 1.04. (Note that a SFE of 0.4 would be a high value for 
an op en cluster, but not necessarily for a UCD, cf. iMurravl 
1200a . Note also that we ignore the probably significant loss 
of stars from the UCD due to the unbinding effect from gas 
expulsion. We return to this in a follow-up paper.) 

The behaviour of a stellar system that loses a large frac- 
tion of its initial mass very much depends on the rate of the 
mass-loss. This behaviour can be characterised by two lim- 
iting cases: 

• Rapid mass-loss (i.e. the mass-loss takes place on a 
timescale shorter or compar able t o the c rossing time): Argu- 
ing with the virial theorem, Irlillsl (|l980l ) finds from analytic 
estimates that a star cluster dissolves if it loses more than 
50% of its initial mass instantan eousl y. Using Af-body in- 
tegrat ions, iBoilv fc Kroupal (|2003r i and lFellhauer fc Kroupal 
l|2005t ) show that the survival of a star cluster is also de- 
pendent on the density profile of the star cluster and the 
velocity distribution of its stars, but a sudden loss of more 
than 67% of its initial mass is critical in any case. 

• Adiabatic mass-loss (i.e. the mass-loss is slow enough 
for the stellar system to stay near virial equilibrium at all 
times): Adiabatic mass- loss does not unbind the remainder 
of the star cluster, but inflates it. The change in radius is 



7-final 
Hnit 



Mjnit 

M flnal 



(16) 



where rinit and Mi n j t are the radius and mass, respectively, 
of the stellar system at the beginning of mass-loss and rfi na i 
and Mfinai the according parameters at the end of mass-loss 
|Kroupall2008T ). 

Important numbers for deciding whether the mass-loss 
from UCDs is rapid or adiabatic are their crossing times, 



which is defined as t c 



2rh/u, with rn being the 



3D half-mass radius and a being the 3 D velocity disper- 
sion of the stellar system l|Kroupal 12001) . Furthermore the 
ratios between half-mass radius and tidal radius, rh/rt, de- 
ci de upon survival of the st ellar system. The data p ublished 
in lEvstigneeva et all l|2007j ) and lHilker et all (120071 ) for the 
properties of UCDs as they are observed today imply rh/rt 
well below 0.1 for most of them (sometimes as low as 0.01) 
and i C r of the order of 1 Myr. This suggests a timespan of 
the order of 40 t cr for the timescale for SN-driven mass-loss. 

In the grid of A^-bod y simulations performed by 
iBaumgardt fc Kroupal (|2007t ). star-clusters are predicted to 
dissolve if they loose 95% of their initial mass, even for the 
most moderate tidal fields (r n /r t = 0.01) and longest du- 
ration for mass-loss (10 t CI ) they consider. However, if the 
stellar system looses only 90% of its initial mass on that 
timescale, the stellar system may not completely dissolve as 
long as the tidal field is weak. The mass fraction of the rem- 
nant of the star cluster that continues to be gravitationally 
bound after it has returned to virial equilibrium is then 0.65 
for rh/rt = 0.01 and 0.35 for r h /r t = 0.033, while the half- 
mass radius increases to approximately ten times its initial 
value in both cases. This implies that the stellar density 
in those systems decreases to less than 10 -3 times its initial 
value. If also the gas leaving intermediate mass stars as they 
evolve into WDs is driven out of the UCDs (e.g. through 
type I SNe), the mass of the UCDs is decreased further to 
0.57 of its original value for 03 = 1.57 and to 0.40 of its 
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original value for Q3 = 1.04. Since this mass-loss would be 
adiabatic, the according change in radius can be calculated 
using eq. (|16[) and is a factor of 1.75 for as — 1.57 and a fac- 
tor of 2.48 for «3 = 1.04. The density would thus be further 
decreased by a factor < 10. At present, the UCDs typically 
have mean central densities from 10 2 to 10 3 Mq pc -3 (fig. 4 
in DHK). A mass-loss of 90% of the initial mass (stellar and 
gas) over the first « 40 t CI and subsequent adiabatic mass- 
loss through the evolution of intermediate-mass stars would 
therefore suggest initial central densities of at least 10 6 to 
10 7 M pc~ 3 . This corresponds to 3.9 x 10 5 to 3.9 x 10 6 stars 
per pc 3 for a 3 = 1.57 and 8.4 x 10 4 to 8.4 x 10 5 stars per 
pc 3 for «3 = 1.04, the fraction of stars more massive than 
8Mq among them being 5.5% for 03 = 1.57 and 23.7% for 
a?3 = 1.04. Such systems would thus have had extensions 
similar to GCs (i.e. r\ of a few pc). Typical total initial stel- 
lar masses would be some 10 7 Mq , implying a population 
of ~ 10 6 stars with m ^ 8Mq for a-j, = 1.57 as well as for 
a 3 = 1.04. 

These numbers underline the extreme nature of UCDs. 
Their stability seems questionable if they would have formed 
with a top-heavy IMF and their contemporary structural pa- 
rameters. However, the smaller extensions and higher masses 
the UCDs must have had before evolutionary processes set 
in imply that the conditions for adiabatic mass-loss were 
fulfilled much better at that time. However, also if mass- 
loss from a UCD with a top-heavy IMF was adiabatic at all 
times, and its stability was therefore not threatened, eq. [15] 
still implies an enormous inflation nd decrease of density for 
it. 

The observation of UCDs today therefore does not con- 
tradict a formation scenario with a very top-heavy IMF for 
them. A more detailed, numerical study of this issue will be 
provided in a follow-up paper. 

5.2 The Star formation rate in UCDs at their 
formation 

The notion that UCDs might be the most massive star 
clusters implies that they formed from a collapsing molec- 
ular cloud. Star formation within the cloud is thought to 
set in as so on as a certain density is reached, which is 
according to iKawamura et all l|l998h at a column density 
m excess of 1.6 x 10 21 N{H 2 ) cm . Defining the size, R, 
of a cloud as (S/71-) ' 5 where S is the total cloud surface 
area, gives typical sizes of 3 p c for the clouds in the sam- 
ple of IKawamura et alj |l998). This corresponds to a mean 
density, p, of « 4Mg pc -3 as the criterion for the onset 
of star formation. It then proceeds rapidly and is com- 
pleted with i n a t i mescale of the order o f a free fall time 
(|Elmegreenl |2000| ; lHartmann et al.l 120011 ). For spherically 
symmetric matter distributions, the free fall time, tg, is 
given as tg = (37r/32Gp) ' 5 , where G is the gravitational 
constant. Note the independence of ts on the total mass. 
The assumption of spherical symmetry for star-forming gas 
clouds thus leads to a time scale of 4Myr on which star 
formation takes place, whereby the bulk of the stars may 
form on an even shorter time scale (f or instance, 80% wit hin 
lMyr in the Orion Nebula Cluster. IProsser et af]|l994h . If 
applied to the UCDs, this suggests that their stellar popu- 
lations formed with star formation rates of « lOMgyr" 1 
to ~ lOOMgyr -1 , depending on the mass of the UCD. 



Given that most if not all stars are formed in star clusters 
IlLada fc Ladall2003h and that the time scale for star clus- 
ter formation appears to be independent of the mass of the 
cluster, these star formation rates would be the highest ever 
to be found in a single star formation event. 



6 SUMMARY AND CONCLUSIONS 

It was shown in previous papers that the dynamical MjLy 
ratios, Ty, of compact stellar systems more massive than 
2 x 10 6 M are not consistent with the predictions from 
simple stellar population models, if the canonical IMF is 
assumed for star formation in them. Out of the possible 
explanations for this result (top-heavy IMF, bottom-heavy 
IMF, dark matter, inaccuracy of the SSP models), the no- 
tion of a top-heavy initial stellar mass function (IMF) in 
dense star-forming regions seems especially attractive. 

With this motivation, we quantify by how much the 
IMF in a sample of massive compact stellar systems (re- 
ferred to as UCDs) has to deviate from the canonical IMF 
in the intermediate and high mass part for the modelled 
Ty to agree with the observed ones. The model constructed 
for this accounts for the different metallicities of the UCDs. 
Several combinations of assumptions concerning age (7Gyr 
or 13 Gyr) and the amount of processed material with very 
high Ty retained by the UCDs besides white dwarfs (no 
remnants of massive stars; 20% of the remnants of massive 
stars; all remnants of massive stars; all material processed 
by burnt-out stars) are considered. The IMF of the UCDs is 
taken to be a three-part power-law that equals the canonical 
IMF below an initial stellar mass of 1 Mq . The exact upper 
mass limit of the IMF (m max = 100 Mq or m max = 150 Mq) 
turns out to have a negligible impact on the results. 

Assuming that all UCDs have the same normalised 
M/L ratio (which is justifiable considering the uncertainties 
of their Ty) and that the processed material retained by the 
UCDs are all white dwarfs and 20% of the remnants of mas- 
sive stars, our model suggests a high- mass IMF slope, a 3 , of 
« 1.6 if the UCDs are 13 Gyr old (i.e. almost as old as the 
Universe) or w 1.0 if the UCDs are 7 Gyr old. If the UCDs 
were assumed to have formed with the canonical IMF, their 
Ty would only be explainable if they contain significant 
amounts of non-baryonic dark matter or dense interstellar 
gas. Note that there would need to be some mechanism that 
inhibits on-going star formation in this case. 

The Ty of the UCDs in the Fornax cluster tend to be 
lower than the ones of the other UCDs. If the Fornax UCDs 
are assumed to be 13 Gyr old they have a normal Ty and 
consequently not at a top-heavy IMF. Assuming that the 
discrepancy between the prediction of SSP models with the 
canonical IMF and the observed M/Lv ratios of the Fornax 
UCDs as high as for the Virgo UCDs suggests an ages around 

7 Gyr for the Fornax UCDs if the Virgo UCDs are taken to 
be 13 Gyr old l|Mieske et alj|2008h . 

The dependency of as on the normalised M/Ly ratio, 
Ty ;I1 , established in eq. (|1 1 p . can be translated into a de- 
pendency of Q3 on the mass of the stellar system, M. This 
is done using the increase of Ty jn with M formulated in 
eq. [T3] and shown in Fig. [3] A possible parametrisation of 
this dependency is given in eq. (|14[) and plotted in Fig. [5] 

The mass-loss due to the evolution of massive stars may 
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reach 90% of the initial stellar mass of a star cluster for very 
top-heavy IMFs, even if primordial gas expulsion is not con- 
sidered. The survival of the UCDs seems not to be threat- 
ened by the mass loss implied by the evolution of the stars 
alone, as this mass loss would be adiabatic. However, the ra- 
diation and evolution of massive stars also drives the expul- 
sion of the primordial gas. The timescale on which this pro- 
cess takes place is critical for the survival of the UCD. This 
is an issue deserving further study. In any case, the results in 
Sections 15.11 and 15.21 suggest that UCDs formed with likely 
central stellar densities of 10 6 to 10 7 Mq pc -3 and possi- 
ble star formation rates of « 10 Mq yr" 1 to fa 100 Mq yr _1 . 
These are among the most extreme sites of star formation. 
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APPENDIX A: THE TOTAL MASS OF THE 
REMNANTS 

We explicitly note the terms with m > m to that arise from 
the integration of the right hand side of eq. ©, if eq. (|HJ is 
inserted for m rcm (m). 

The contribution of the white dwarfs to M m , M m ,wD, 
can be written as 

-Mm,WD,v _ 0.109 ks , 2-o 3 2-OJ3N 

M 2-aa 

+ 0394fc3 x(8 i-3_ i- 3)i 

1 — Q 3 

with the masses of the stars in Solar units, mto was argued 
to be « 1 M in Section 13X21 

The contribution of neutron stars to M m , M mi NS, can 
be written as 

A^ns = L35fc3 x (25l _« 3 _ gl - a3) _ (A2) 

Mr 1 - 03 



